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ABSTRACT 

There still exists the possibility that the phenomenon of gamma-ray bursts (GRBs) is a mixture 
of events of different nature, even within the class of long (>2 s) bursts. We try to put statistical 
constraints on a possible non-cosmological component using the uniform GRB catalog of 
Stern &Tikhomirova* obtained from an overall scan of the full 9. 1 year BATSE 1024 ms data. 
The sample consists of 3906 GRBs and includes non-triggered bursts with peak fluxes down 
to 0. 1 photons cm^^ s^^ . We find no significant deviations from isotropy. The constraints on a 
non-cosmological population are still weak. The allowed contribution of a GRB subpopulation 
originating from an extended galactic halo is ^ 60% and the upper limit on an Eucledian 
component (e.g., nearby galactic or non-cosmological extragalactic GRBs) is 23%. The results 
concern mainly the class of long GRBs. 
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1 INTRODUCTION 

At present there are direct redshift measurements for 23 Gamma- 
Ray Bursts'^). All of them are cosmological. If we believe that 
all GRBs are of the same nature, then we must reject all non- 
cosmological models of GRBs and to close this issue. 

However, first we should mention that long (>2 s) and short 
(<2 s) GRBs are probably events of different nature (Kouveliotou 
et al. 1993) and all of the redshift measurements have been done 
only for long bursts. Moreover, the impression of the uniformity 
of GRBs could arise from their diversity whereas they can be a 
mixture of events of different origin even within the class of long 
bursts. For example, Horvath (1998) suggested the existance of an 
"intermediate" class of GRBs with durations of several seconds and 
softer spectra. The evidence for this is still not statistically convinc- 
ing. Nevertheless such a possibility should not be discarded. 

In this work, we revisit statistical studies of non-cosmological 
models which were done before the discovery of afterglows. Now 
we have a much larger statistics which will not be extended during 
the next few years. Therefore it is interesting to outline constraints 
on non-cosmological (i.e., local galactic, galactic halo, or low red- 
shift extragalactic) GRB subpopulations which we can obtain from 
the present GRB sample. We concentrate only on those models 
which can provide a reasonable isotropy of GRBs and therefore 
we do not consider the galactic disk population (except for a very 
local one). 



* E-mail: jana@anubis.asc.rssi.ru (YT); stern@lukash.asc.rssi.ru (BES); 
svensson@astr o su.se fRS) 
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Previous attempts to constrain Galactic GRB scenarios (see, 
for example, Hakkila et al., 1994, Loredo & Wasserman, 1998) 
showed that the least constrained local model is an extended galac- 
tic halo of bursters. According to Hakkila et al. (1994), for some 
narrow range of parameters, all GRBs could originate from a galac- 
tic halo being still in agreement with the angular and the bright- 
ness distributions of GRBs as observed by the Burst And Transient 
Source Experiment (BATSE)(Fishman 1992). 

It is also interesting to recalculate a possible Euclidean com- 
ponent (i.e., a population homogeneously distributed in Euclidean 
space). It could be represented by a local Galactic population at dis- 
tances <500 pc or by extragalactic sources at z^I. The estimate 
of Kommers et al. (2000) for the BATSE sample with peak fluxes 
down to ~ 0.2 photons cm~^ s~^ is about 10%. 

The deepest and largest GRB sample was found by search for 
non-triggered GRBs in the BATSE continuous records of Stern et 
al. (2000, 2001). Their new uniform catalog* (UC) includes 3906 
(2068 triggered and 1838 non-triggered) events with peak fluxes 
down to ~ 0. 1 photons cm~^ s~^ (a factor 2 lower than the BATSE 
trigger threshold). Moreover, they applied a new method to measure 
the efficiency of the search which is important for the reconstruc- 
tion of the real log A^-log P distribution. 

Note, that all bursts of the UC were extracted from the DIS- 
CLA BATSE data of the 1024 ms time resolution. The bursts which 
can be classified as short (I bin events) are about 12% in the sam- 
ple. Therefore, the constraints will concern mainly the class of 
long GRBs. In §2, we present the results of a general test of the UC 
for large scale isotropy. We constrain the possible contribution of a 
halo population in §3, and of an Euclidean component in the GRB 
statistics in S4. 
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Figure 1. Sky distribution of the 3906 GRBs in tlie UC on an Aitoff- 
Hammer projection in Galactic coordinates. 



2 TEST FOR ISOTROPY 

The preliminary test of the UC for isotropy was done in 
Tikhomirova & Stern (2000). Here, we extend the analysis to the 
complete sample. 



2.1 Statistical tests 

Figure 1 shows the distribution of the locations of all UC GRBs in 
galactic coordinates. We use the following tests for the large scale 
isotropy (see Hartmann & Epstein 1989 for test 1) and Briggs 1993, 
1994 for tests 2)-5)) : 
1) The dipole vector 



7? 



iV 



where rt is the unit vector in the direction of a GRB and TV is the 
number of bursts in the sample. ^ characterizes the dipole moment 
in a coordinate free way and tests for large scale anisotropy. 

2) (cos(6)), where Q is the angle between the Galactic center and a 
GRB. This quantity characterizes the dipole moment in the Galac- 
tic system and tests for a Galactic population. 

3) {sin'^(6) — 1/3), where h is the Galactic latitute of a burst. This 
quantity characterizes the quadrupole moment in the Galactic sys- 
tem and tests for a local Galactic population. 

4) {sin((5)), where & is the declination of a burst, and 

5) {sin^(5) — 1/3), which characterizes the dipole and quadrupole 
moments in the equatorial system and tests for systematic errors. 

The statistical errors are a — (l/\/3iV) for the components 
of ~k and the tests 2) and 4); and a = (\/4/V45Af) for the tests 
3) and 5). The statistical error for the absolute value of ^ is cr = 
(l/\/7V). The total error has an additional component associated 
with the errors of the locations of individual GRBs. Although the 
location accuracy for the UC GRBs is as low as several degrees 
(Tikhomirova & Stem 2000) this component is smaller by an order 
of magnitude than the statistical one as was checked using Monte- 
Carlo simulations. 



2.2 Expected values 

The values of j ^ | and of the statistics 2)-5) are zero for an isotropic 
distribution of bursts observed by an ideal instrument. The values 
corrected for the non-uniform BATSE sky exposure (Paciesas et al. 
1999) are listed in the second column of Table 1. In the BATSE 



catalogs, these quantities usually have their expected values (see 
Paciesas et al. 1999). However, other systematic effects may ap- 
pear as a result of a non-uniform burst selection. One such effect is 
caused by an active CygX- 1 , which reduces the efficiency of weak 
GRBs detected in that direction. 

Such effects may be revealed by the test bursts used in the 
scan of Stern et al. (2000, 2001) to estimate the scan efficiency. 
Test bursts are initially isotropic and are subjected to the same sys- 
tematic effects as real bursts. However, their statistics is limited 
(~ 5300). Therefore, we consider both the values expected for the 
BATSE sky exposure (column 2 in Table 1) and the values observed 
for test bursts (column 4 in Table 1). 



2.3 Results 

The results of the tests together with the la statistical errors are 
listed in columns 4-7 of Table 1 . The values for the 4Br catalog are 
given for comparison in column 3. (Note, that ^ was not given in 
the 4Br catalog.) Test bursts show a marginally significant dipole 
(at the ~ 2.5(7 level) in the direction opposite to Cyg X-1, which is 
caused by a strong deficit of detected weak bursts in a cone ~ 25° 
around Cyg X-1. A similar effect appears for the sample of real 
GRBs. 

Table 2 shows the deviations from the values expected for the 
BATSE exposure (columns 2-5) and from the values for test bursts 
(columns 6-8). The deviations from the value for the test bursts are 
given in units of a determined as 



real bursts ' 



(1) 



All deviations are within 1.5cr so the results of the tests are consis- 
tent with isotropy. 

We also checked the excess of bursts towards M31 for the full 
sample in the UC as well as for only weak (<0.4 photons cm~^ 
s~^) bursts. There is no excess. 



3 TEST FOR AN EXTENDED GALACTIC HALO 
SUBPOPULATION OF GRB SOURCES 

3.1 The Model for an Extended Halo of GRB sources 

The traditional model of a halo is associated with old neutron 
stars which have been ejected from the galactic disk (Shklovskii 
& Mitrofanov 1985). The test for an extended halo population is 
based on the displacement of the Solar system from the Galactic 
center (8.5 kpc) and the presence of the nearby galaxy M31 (670 
kpc away) which should have a similar extended halo. Then we 
must check two things: the dipole moment in the direction of the 
Galactic center and the excess of bursts around M31. The larger the 
halo of bursters, the smaller is the dipole moment but then GRBs 
from M3 1 become more visible. The same logic was used in previ- 
ous works (e.g., Hakkila et al., 1994). 

In the approximation of a steady outflow of neutron stars, we 
have a spatial distribution of bursters with a Coulomb-like tail and 
a core: 



Nir) = Norl/ir^,+r^), 



(2) 



where N{r) is the burster density per unit volume, A'o is the density 
at the center, r is the distance from the Galactic center, and Vc is 
the core size. The distribution should decline faster than 1/r^ at 
some r because a neutron star cannot emit GRBs for an infinite 
time or just because of the finite age of the Galaxy. Expression (0) 
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Table 1. Results of the large scale isotropy test. 





expected 


4Br 






UC 




statistics 


for 


catalog 


test bursts 


all bursts 


triggered 


non-triggered 




isotropy 


(1637)* 


(5345)* 


(3906)* 


(2068)* 


(1838)* 


17?1 


0.018 


0.029 


0.026 


0.021 


0.020 


0.030 






±0.025 


±0.014 


±0.016 


±0.022 


±0.023 


Rl 


0.000 


0.005 


-0.022 


-0.012 


-0.006 


-0.019 






±0.014 


±0.008 


±0.009 


±0.013 


±0.013 


R2 


0.000 


0.015 


0.014 


0.014 


0.005 


0.024 






±0.014 


±0.008 


±0.009 


±0.013 


±0.013 


R3 


0.018 


0.024 


0.003 


0.010 


0.019 


-0.000 






±0.014 


±0.008 


±0.009 


±0.013 


±0.013 


(cos(e)> 


-0.009 


-0.025 


-0.012 


-0.016 


-0.013 


-0.019 






±0.014 


±0.008 


±0.009 


±0.013 


±0.013 


(sin2(6) - 1/3) 


-0.004 


-0.001 


0.001 


-0.007 


-0.007 


-0.006 






±0.007 


±0.004 


±0.005 


±0.007 


±0.007 


(sin(5)> 


0.018 


0.024 


0.003 


0.010 


0.019 


-0.000 






±0.014 


±0.008 


±0.009 


±0.013 


±0.013 


(sin2(5) - 1/3) 


0.024 


0.025 


0.022 


0.025 


0.024 


0.026 






±0.007 


±0.004 


±0.005 


±0.007 


±0.007 



* the number of events in each sample 



Table 2. Deviations from the expected values for large scale isotropy. 



deviation from the value deviation from the value 

statistics expected for the isotropy for the test bursts 





4Br 




UC 






UC 






catalog 


all 


trig. 


non-tr 


all 


trig. 


non-tr 


it* 


0.7(7 


1.2cr 


0.3(7 


1.5(7 


0.6(7 


0.9(7 


0.4(7 


(cos(e)) 


-l.lcr** 


-0.7a 


-0.3(7 


-0.8(7 


-0.3(7 


-0.1(7 


-0.5(7 


{sin2(6) -1/3) 


±0.4cr** 


-0.5a 


-0.5(7 


-0.2(7 


-1.1(7 


-1.0(7 


-0.8(7 


(sin(5)) 


+0.4(7** 


-0.9a 


+ 0.1(7 


-1.3(7 


+0.6(7 


+ 1.1(7 


-0.2(7 


(sin2(5) - 1/3) 


+0.1(T** 


+0.2a 


+ 0.0(7 


+0.3(7 


+0.5(7 


+0.3(7 


+0.5(7 



* deviation of the vector 



** Paciesas et al. (1999) 



gives a suitable asymptotic slope of the log A*' - log P distribution: 
dN/d{\ogP) oc P~°' with a = 0.5 for small brightnesses P in 
agreement with the observed data (Kommers et al. 2000, Stern et 
al. 2001). 

We assumed identical halos for both our Galaxy and M31 and 
tried two cutoffs for the distribution (^: at 300 kpc (i.e., halfway 
between the Galaxy and M31) and at 800 kpc which could rep- 
resent the case of an infinite distribution (we are within the M31 
halo in this case). This rough model seems to be sufficient for an 
approximate estimate. 

The GRB luminosity function was assumed to be a lognormal 
distribution: 



dN/dL = exp[- log\L/Lo)/al], (3) 

where L is the GRB peak luminosity, Lo is the average GRB peak 
luminosity, and ctl is the width of the GRB luminosity function. 

The model distributions of GRBs were obtained using Monte- 
Carlo simulations. The variable parameters were: the core radius, 
Vc, the average GRB peak luminosity, Lo, and the width of the GRB 
luminosity function, ctl- The averaged peak luminosity Lo was ex- 
pressed as Pioo measured in units of 2- 10*^ erg s~^ in the 50 - 300 
keV energy range. This is the luminosity which produces a flux of 
1 photon s"'^ cm~^ at a distance of 100 kpc (Lo ^ 2 x iC'^Pioo). 
The spatial distribution of GRBs was sampled according to 
, where r is the distance of a GRB from the center of the Galaxy 
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or from the center of M31 (with equal probabihties). The locations 
of GRBs on the sky were sampled imposing experimental loca- 
tion errors depending on the "observed" brightness. The intrinsic 
peak brightness of a GRB was sampled according to equation (^. 
The "observed" distribution of simulated GRBs was folded with the 
BATSE exposure function (Paciesas et al., 1999) and with the de- 
tection efficiency determined using the test bursts method of Stem 
etal. (2000, 2001). 



3.2 The constraints on the halo subpopulation 

We represented the full sample of GRBs as — l)So +^Sh where 
5*0 is an isotropic subsample (presumably cosmological) with an 
unknown log - log P distribution, and Sh is the model halo 
subsample. We want to determine the value of the halo fraction ^ 
for which the following four criteria are satisfied: 

(i) The dipole moment, Dg, towards the Galactic center should 
be within the observed 2a upper limit: Cg — 0.009 

(ii) The fraction of bursts with locations within 25° from the 
M31 location, A^25, should not exceed the observed 2a upper limit 

which is C25 = (A^rcal - iVoxp + 2cr)/iVrcal = 0.14 

(iii) The same should apply for an 18° region around M31, for 
Ari8:Ci8 =0.21 

(iv) The log A'^ - log P distribution should be consistent with the 
observed one: 

- the sum of the distributions for subsamples So and Sh should give 
the observed one, 

- the distribution of the subsample So should be smooth, i.e., the 
difference between neighbouring points should not exceed 3a, and 

- the distribution of So should not bend down too sharply at low 
brightnesses, i.e., not steeper than dN/d{\og P) oc p+"'5 (see Fig. 
2). This limiting slope corresponds to a reasonably sharp turnover 
of the log A'^-log P distribution at its dim end. 

For the observed distribution, we used the log A'^ - log P curve 
corrected for the detection efficiency of the scan of Stem et al. 
(2001). All limits were calculated for the full sample of the UC 
(3906 GRBs with peak fluxes down to 0.1 photons cm~^ s~^). 

We searched for the fraction ^ at which the full model sample 
of GRBs satisfies all four criteria. For the first three criteria, ^ is 
just the minimum of the ratios Cg / Dg , C25 /N25 , and Cis /Nis . 

We also tried to apply criteria 2 and 3 for weak bursts only 
(with peak fluxes <0.4 photons cm^^ s^^). This case does not give 
stronger constraints. 



3.3 Results 

Comparing results for different halo cutoffs, we found that all con- 
straints are slightly stronger for the 300 kpc cutoff. Here we present 
the results for the 800 kpc cutoff only, as it represents the case of 
an infinite halo and gives more conservative constraints. 

The hypothesis that all GRBs originate from a Galactic halo is 
inconsistent with the data for any parameters. However, this fact 
is only of academic interest as observations of GRB afterglows 
showed that at least a substantial part of GRBs has a cosmologi- 
cal origin. 

The results of the Monte-Carlo simulations for different pa- 
rameters: Vc, see equation (^, Pioo, and ctl, see equation (^, are 
shown in Figure 2 as maps of isocontours of the allowed halo frac- 
tion 5. The highest allowed fraction ^ is 0.63 for aL = 1 and 0.55 
for aL ~ 2 (for ctl = 0, ^ = 0.7, but the case of ctl = 0, i.e.. 
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Figure 2. Maps of constant levels of the allowed fraction ^ of the halo 
subpopulation in the whole sample of the UC: a) for ctl = 1, b) for (Tl = 2. 

is the halo core radius, and Pioo is the GRB peak count rate at a distance 
of 100 kpc. 



for a standard candle GRB luminosity, is hardly realistic, and we 
do not discuss it further). The brightness intervals where the frac- 
tion of the galactic subpopulation can exceed 0.5 are limited to: 
0.3 < Pioo < 1-0 for CTL = 1 and 0.012 < Pioo < 0.03 for 
(jl = 2. In terms of the average peak luminosity these intervals 
are (1.0 - 3.3) X 10*^ ergs"^ and (1.8 - 4.5)x 10"*° erg s"\ respec- 
tively (note the difference between arithmetical mean and logarith- 
mic mean for the case if lognormal distribution). 

The total rate of GRB events in the halo required to reproduce 
the observed GRB data depends on ctl. In the case of ctl = 1, 
the rate of GRBs out to 800 kpc should be ~ 8000 per year 
(TVo = 0.004 kpc"^). In the case (Tl = 2, this number is ~ 45 
000 (TVo = 0.02 kpc~^). Nevertheless, the total power emitted in 
the form of GRBs implied in the two cases is almost the same: 
erg year , where T is the average duration of GRBs 
(a reasonable estimate is T ~ 10 s). In both cases, M31 is at the 
edge of the sampling volume. The concentration of GRBs around 
the M31 location should be still visible for ah ~ 2 and hardly 
visible for ctl = 1 (see Fig. 3). 



4 CONSTRAINTS ON THE EUCLIDEAN 
SUBPOPULATION 

We used the same logic as in §3 to constrain the possible fraction 
of a GRB subpopulation homogeneously distributed in space (the 
Euclidean component). 

We represent the whole sample of GRBs as — 1)5*0 + ^Se, 
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Figure 3. The simulated distributions for the case of maximum halo fraction ^: a) for ctl = 1, b) for ctl = 2. Panels on the left, the log N - log P 
distributions: the empty squares represent the observed data according to the UC, the filled squares represent the halo population, and the crosses correspond 
to the cosmological population. Panel on the right, the sky distribution: the whole sample (cosmological + halo). The halo of M31 is visible only for (Tl = 2. 
M31 is located at I « 122° and 6 » —21°. on an Aitoff-Hammer projection in Galactic coordinates 



where 5*0 is an isotropic subsample (presumably cosmological) 
with an unknown log N - log P distribution, Se is the homoge- 
neous (in Euclidean space) subsample, and ^ is the possible frac- 
tion of GRBs of the homogeneous subsample in the full sample. 
For constraining ^, we used only criterion 4 from §3.2. 

The upper limit on an Euclidean component is 23% for the 
full sample of the UC (3906 GRBs with peak fluxes down to 0.1 
photons cm~^ s~^). For the sample from the UC with peak fluxes 
down to the BATSE trigger threshold (« 0.2 photons cm~^ s~^), 
the limit is about 14%. The same estimate of Kommcrs et al. (2000) 
is about 10%. They used their log - log P distribution with fluxes 
down to about 0.18-0.20 photons cm~^ s~^ and the difference in 
the constraints on the Euclidean component results from the dif- 
ference in the estimates of log A'^ - log P distribution. Kommers' 
estimate for the number of dim GRBs is lower. 

This constraint on the Euclidean component should be applied 
to any kind of local Galactic GRB subpopulation (r < 300 pc) and 
to extragalactic sources at small redshifts {z <^ 1). If the sample 
of GRBs has a ~ 20% fraction of an Euclidean subpopulation, then 
this population dominates among weak GRBs. If this subpopulation 
is intrinsically different, it can be revealed in correlations between 
some GRB properties and the brightness. 



5 CONCLUSIONS 

The general results of our analysis of the new uniform catalog of 
GRBs can be formulated as follows: 

-the large scale isotropy of the GRB sky distribution is con- 



firmed for the larger and deeper sample of 3906 BATSE GRBs with 
peak fluxes down to ~ 0.1 photons cm~^ s~^; 

-the possible fraction of GRBs from an extended galactic halo 
can be up to ~60%. This is, however, only for narrow intervals of 
the parameters; 

-the Euclidean component limited by a fraction of GRBs 
above the threshold 0.1 photons cm~^ s~^ is about 23%. 

These constraints are still weak, so a considerable fraction 
of even long GRBs can be prescribed to some kind of a non- 
cosmological source population. Stronger constraints can hardly be 
obtained within the framework of the statistical approach. We be- 
Ueve that tighter constraints can be imposed from future afterglow 
observations. 
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